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1. Introduction 

Exactly solvable examples played an important role in the understanding of VT 
symmetric quantum mechanics and its unusual features. The first VT symmetric 
potentials, i.e. potentials which are invariant under the simultaneous action of the 
V space- and T time reflection operations were found numerically The most 

surprising result was that these one-dimensional complex potentials possessed real 
energy eigenvalues, which however, turned pairwise into complex conjugated pairs 
as some potential parameter was tuned. This mechanism was interpreted as the 
spontaneous breakdown of VT symmetry, since the potential remained VT invariant 
throughout, but the eigenf unctions associated with the discrete spectrum lost this 
property as the spectrum turned into complex gradually with the tuning of the potential 
parameter. Further examples have been found in numerical j2], semiclassical |3j and 
perturbative [I] studies, but a number of quasi-exactly solvable (QES) [5] and exactly 
solvable [01 IB] VT symmetric potentials have also been identified. These latter 
potentials were analogues of Hermitian exactly solvable potentials with the property 
that their real and imaginary components were even and odd functions of the coordinate 
x, respectively. 

More recently VT symmetric problems have been analysed in terms of pseudo- 
Hermiticity [H], and their unusual features have been interpreted in terms of this more 
general context. A Hamiltonian is said to be ?7-pseudo-Hermitian if there exists a linear, 
invertible, Hermitian operator for which W = nHn" 1 holds. It has been shown that 
for systems with Hamiltonians of the type H = p 2 + V(x) VT symmetry is equivalent 
with P-pseudo-Hermiticity. With other choices of rj further complexified Hamiltonians 
can be generated, which might also have real energy eigenvalues but do not fulfill 
VT symmetry [TUJ [n]. Furthermore, with n = 1 ^-pseudo-Hermiticity reduces to 
conventional Hermiticity. 

Here we restrict our analysis to VT symmetric problems, and in particular, to 
exactly solvable ones. A number of peculiar features of VT symmetric potentials became 
apparent only during the analysis of these potentials. 

• It has been observed that several exactly solvable VT symmetric potentials possess 
two sets of normalisable solutions [01 EHl E2j in the sense that there can be two 
normalisable states with the same principal quantum number n. The second set 
of solutions can appear in several ways. For potentials which are singular at 
the origin the problem can be redefined on various trajectories of the complex 
x plane such that the integration path avoids the origin and the solutions remain 
asymptotically normalisable. (This latter feature is similar to some numerically 
solvable VT symmetric problems pQ.) A special case of this scenario is obtained 
when an imaginary coordinate shift x — > x + is is employed [7] . The advantage of 
this scenario is that the discussion of the VT symmetric potential remains rather 
similar to that of its Hermitian counterpart, and by suitable and straightforward 
modification of the formalism it can also be interpreted as a conventional complex 
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potential denned on the x axis. This imaginary coordinate shift can be employed 
[7j to all the shape-invariant potentials (e.g. the VT symmetric harmonic 
oscillator [Oj or the generalised Poschl-Teller potential ^1]) with the exception of 
the VT symmetric Morse jT3] and Coulomb JH] potentials. The cancellation of the 
singularity then regularizes the solution which would be irregular at the origin in 
the Hermitian setting. A different mechanism appears for potentials which are not 
singular in their Hermitian version, such as the VT symmetric Scarf II potential, 
which is defined on the full x axis. In this case the second set of normalisable 
solutions originates from states which have complex eigen-energy in the Hermitian 
case, but which turn into normalisable states with real energy when the potential is 
forced to become VT symmetric and the VT symmetry is not broken spontaneously 
[12] H3 [TJJ HH]. The two set of solutions are distinguished by the quasi-parity 
quantum number [TTij . 

• In the process of generating the spontaneous breakdown of VT symmetry by 
tuning the potential parameters it was found that the pairwise merging of the 
energy eigevalues and their re-emergence as complex conjugated pairs occurs at the 
same value of the potential parameter HZ] • In other words, the spontaneous 
breakdown of VT symmetry is realised suddenly in the case of shape-invariant 
potentials, as opposed to a gradual process observed in the case of numerical 
examples [T] l2*T]. 

• For certain exactly solvable (and shape-invariant) examples, such as the VT 
symmetric harmonic oscillator [22J and the Scarf II potential j2B] it was found 
that there are two 'fermionic' SUSY partners of the original 'bosonic' potential, and 
they are distinguished by the quasi-parity quantum number carried by the 'bosonic' 
bound states. (In the latter case this has also been found in other realisations of 
SUSYQM •) This doubling of the partner potentials is an obvious consequence 
of the fact that there are two nodeless normalisable solutions corresponding to the 
'ground state' in the two segments of the spectrum with quasi-parity q — +1 and 
— 1. Furthermore, it was also established that the 'fermionic' partner potentials 
are VT symmetric themselves too, in case the 'bosonic' potential has unbroken 
VT symmetry, while they cease to be VT symmetric if the VT symmetry of the 
'bosonic' potential is spontaneously broken. 

The peculiar features mentioned above have been observed until now for the VT 
symmetric version of shape-invariant potentials, while potentials beyond this class 
often behaved in a different way. This naturally raises the question whether these 
features also characterise non-shape-invariant, but exactly solvable examples. Natural 
candidates for this analysis are Natanzon-class potentials [2H] which have the property 
that their bound-state solutions are written in terms of a single hypergeometric or 
confluent hypergeometric function. This potential class depends on six parameters, but 
it is prohibitively complicated in its general form, so its subclasses with two to four 
parameters have been analysed in detail until now j2El EH EE1 EE1 EDI EI] • 
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Perhaps the most well-known member of the Natanzon-class is the Ginocchio 
potential, which has a one-dimesional version defined on the x axis and a radial 
one, which has an r _2 -like singularity at the origin j2Z|. An important feature of 
this potential is that for a special choice of a potential parameter it reduces to a 
shape-invariant potential, namely to the Poschl-Teller hole (in one dimension) and 
to the generalised Poschl-Teller potential (in the radial case). The one-dimensional 
version of the Ginocchio potential has been analysed in an algebraic framework, and an 
su(l,l) algebra has been associated with it 32., the discrete non-unitary irreducible 
representations of which correspond to resonances in the transmission coefficients. 
(These states have been identified as quasi-bound states in an independent study [S3]-) 
Furthermore, it was also shown that this algebra reduces in two different shape-invariant 
limits to an su(l,l) potential algebra and to an su(2) spectrum generating algebra [34J. 
The phase-equivalent supersymmetric partners of the generalised Ginocchio potential 
have also been derived in a completely analytic form |35j . 

Although the generalised Ginocchio potential is an 'implicit' potential, i.e. the 
z(r) function which is used to transform the Schrodinger equation into the differential 
equation of the hypergeometric function F(a, b; c; z) is known only in an implicit form as 
r(z), nevertheless, similarly to a number of other 'implicit' potentials |25] . this does not 
restrict the applicability of the formulae, because V(r) and the wavefunctions can be 
determined to any desired accuracy, and all the calculations involving these quantities 
(matrix elements, etc.) can be evaluated analytically [33]. Furthermore, we shall see 
that the imaginary coordinate shift which is essential to impose VT symmetry on the 
generalised Ginocchio potential can also be implemented without complications. 

We note that a Natanzon-class potential, the so-called DKV potential has already 
been analysed in the VT symmetric setting by the point canonical transformation 
of a shape-invariant potential [3^], but it was found that it has to be defined on a 
curved integration path. Nevertheless, it also showed similarities with VT symmetric 
shape-invariant potentials, as its spectrum was also richer that that of its Hermitian 
counterpart. 

In section [2 we present the Hermitian version of the generalised Ginocchio potential 
for reference, and in section El we construct its VT symmetric version. Section H] deals 
with the supersymmetric aspects of this potential, while in sectional a summary of the 
results is presented. 

2. The generalised Ginocchio potential 

The first version of the Ginocchio potential was introduced as a one- dimensional 
quantum mechanical problem which is symmetric with respect to the x — > — x 
transformation [26 . Later it was generalised to a radial problem j2Z]> which also contains 
an r _2 -like singular term at the origin. (This latter version of the potential also allows 
a particular functional form of an effective mass, but it can be reduced to a constant 
value by setting one of the parameters (a) to zero.) Following the notation of Ref. [SS] 
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we define the generalised Ginocchio potential as 



V{r) 



7 4 (g(g + l) + l-7 2 

7 2 + sinh 2 u 
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+ 7 4 A(A- 1) 



+ 
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coth u 
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(1) 



4(7 2 + sinh 2 u) 2 4(7 2 + sinh 2 u) 3 

where we changed the notation of Ref. [23 to make it more suitable for our purposes. 
This form can be obtained from the original formulae by setting a — 0, ai — A — |, 
V[ — s, P n i = fi, A = 7 and y = sinhw^ 2 + sinh 2 u)~^ . 

The (generalised) Ginocchio potential is an example for 'implicit' potentials, 
because it is expressed in terms of a function u(r) which is known only in the implicit 
r(u) form: 

1 



T 



tanh f (7 + sinh u) 2 s i n h u 



+ (7 2 - 1) 2 tan 1 ((7 2 - 1) 2 (7 2 + sinh 2 u) 2 sinh 



u 



(2) 



r can take values from the positive half axis, which is mapped by the monotonously 
increasing implicit u(r) function onto itself. This function is, actually, the solution of 
an ordinary first-order differential equation 
du 7 2 cosh u 
dr 



(3) 



(7 2 + sinh 2 u) 2 

defining a variable transformation connecting the Schrodinger equation with the 
differential equation of the Jacobi (and Gegenbauer) polynomials j2Zj- It can be seen 
from Eqs. (j2J) and Q that u(r) behaves approximately as jr near the origin, and as 
7 2 r for large values of r. In the 7 — > 1 limit u becomes identical with r, and (0) reduces 
to the generalised Poschl-Teller potential. 
Bound states are located at 



4 2 

"7 V n » 



where n varies from to n max defined below and 



Y 



-(2n + A + ^) + [(2n + A + ^) 2 (l- 7 



+7 2 (s+^) 2 ; 



(4) 



(5) 



All the terms in are finite at the origin, with the exception of the last one, which 
shows r _2 -like singularity there, and can be considered either as an approximation of 
the centrifugal term with I — A — 1 (A integer), or as a part of a singular potential with 
arbitrary / 7^ A — 1. Setting A = 1 or we get the 'simple' Ginocchio potential pp] 
defined on the line. 

The bound-state wavefunctions are expressed in terms of Jacobi polynomials 
^ n ( r ) = A4( 7 2 + sinh 2 u)* (sinh u)\coshu)-^- x -^ P^ n ' X '^ ] {2 tanh 2 u- 1) (6) 
which reduce to Gegenbauer polynomials [SZj for A = 1. The normalisation is given by 

2 7 2 n! T(fi n + A + n + + A + 2n + \) 



A4 



r(/i n + n + l)r(A + n + !)0u n 7 2 + A + 2n + \) 



(7) 
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Considering that the r — > oo asymptotical limit corresponds to u — > oo (see Eq. (J2J)), 
the wavefunctions become zero asymptotically if fi n > holds. Applying this condition 
to Eq. (jHJ) we find that the number of bound states is set by n max < |(s — A). 

3. symmetrisation of the generalised Ginocchio potential 

The first step in the VT symmetrisation of the generalised Ginocchio potential is 
performing the imaginary coordinate shift which allows its extension to the full x axis 
by cancelling the singularity at the origin. This imaginary coordinate shift is a constant 
of integration from (jHJ), and it modifies (J2J) such that r — > x + is. Here we also switched 
to x instead of r to indicate that the original radial potential is extended also to the 
negative x axis, following the standard treatment of VT symmetric potentials. Similarly 
to the Hermitian case, the variable transformation is determined by an implicit formula, 

x + is = — tanh -1 {{j 2 + sinh 2 u)" 5 sinh uj 

+ (7 2 - 1)2 tan" 1 ((7 2 - 1) 3 (7 2 + sinh 2 u)~2 sinh uj (8) 

however, now u takes on complex values. Figure ^ shows the u(x) function for a 
particular value of 7 and e. This function varies smoothly and it is odd under the VT 
transformation: VTu(x) = —u(x), i.e. its real and imaginary components are odd and 
even functions of x, respectively. Asymptotically the relation u(x) — > x -^±<x J 2 (x + ie) 
holds, while near x = there is a 'kink' in both the real and the imaginary component 
of u(x). 

The VT transform of sinh u(x+ie) is — sinh u(x+ie) (which can be seen analytically 
too by series expansion), and this also determines the VT transform of the potential 
(0). In particular, we may note that it contains sinh-u everywhere as sinh 2 -u (including 
also the term with coth 2 ), so (0) is VT symmetric if all the coupling coefficients are 
real. This restricts 7 2 , s(s + 1) and A(A — 1) to real values. The latter two requirements 
allow the following values of s and A: 

real when s(s + 1) > — \ ( real when A(A — 1) > — \ 

1 A = < 1 (9) 

— +ia when s(s + 1) < -\ i2 + ^ when A(A - 1) < 

We are going to discuss these possibilities later. 

Let us now analyse the two independent solutions of the generalised Ginocchio 
potential, written in the form hypergeometric functions. From among the possible 
linear combinations we chose the ones which have different behaviour at the origin in 
the e — > limit: 

foix) ~ (7 2 +sinh 2 w) 1/4 (coshM) a - fo (sinhM) a+6 " c+ ^F(a,a-c+l;a+6-c+l;-sinh 2 M) ,(10) 
M x ) ~ (7 2 +sinh 2 ii) 1/4 (coshM) a - b (sinhM) c - a - 6+ 5F(l-fe,c-6; c-a-b+1; - sinh 2 u) , (11) 
where a, b and c have to satisfy the following relations: 



1 ± fi a + b — c = ±(A — -) a — b = ± 



1/2 



±uj .(12) 
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With the conditions (J 12)) the two independent solutions can be written as 
^i(x) ~ (7 2 + sinh 2 w) 1/4 (cosh w) ±aJ (sinh w) A 

x F(i( / u + A + ^±^),i(- y u + A + ^±u;);A + i;-sinh 2 M) , (13) 

■0 2 (x) ~ ( 7 2 + sinh 2 w) 1/4 (coshM) ±a; (sinhn) 1_A 

x ^(-(—fj, - A + ^ ± w), ^(/x - A + ^ ± w); ^ - A; - sinh 2 u) . (14) 

Note that the same two functions are obtained irrespective of the signs chosen in the 
first two equations in ()12j) . while the sign of u> remains to be determined from the 
normalisability conditions of the wavefunctions. 

Up to this point the functions (fT3J) and (fTl^ supply the general solutions for the 
energy eigenvalue E = — 7 4 /U 2 . In order to obtain solutions belonging to discrete energy 
eigenvalues one has to set one of the first two arguments of the hypergeometric functions 
to the non-positive integer value —n, reducing them to Jacobi polynomials |HZj • We find 
that in contrast with the Hermitian case, normalisable solutions can be obtained in two 
different ways, corresponding to the condition 

2n + l+ fi nq + q{\-^) -u = , (15) 

where q = 1 and q — — 1 holds for (|13j) and (|14|) . respectively. In this case the two 
solutions can be written in a compact form as 

^nq{x) ~ (7 2 + sinh 2 M) 1 / 4 (cOshM)- 2n - 1 -^-^ A -5)( s i n h M )5+9(A-i) 

x p ( 9 (A-I),- 2 n-l- M (A-I)) (cosh(2w)) (16) 

Here n is the principal quantum number labelling the bound states and q is the quasi- 
parity q = ±1 JQ. This quantum number characterises the solutions of VT symmetric 
potentials, but the potential itself does not depend on it. Its name originates from the 
analysis of the VT symmetric version of the one-dimensional harmonic oscillator: in 
the Hermitian limit of this potential (i.e. for e —>■ it essentially reduces to the parity 
quantum number. 

We see that the two solutions are distinguished by the q = ±1 quasi-parity quantum 
number, similarly to potentials belonging to the shape-invariant class. Actually, the 
corresponding solutions of the VT symmetric generalised Poschl-Teller potential ^1] 
can be obtained by setting 7 = 1. It is also obvious that normalisability requires 
Re(/z ng ) > 0. 

In order to obtain explicit expression for fi nq one has to combine ()15j) with ()12)) and 
to solve a quadratic algebraic equation for /i = /i ng : 



1 

^1 ~ ^2 



2 n + i + g (A-i)) + 7 2 ( s + i) 2 + (i- 7 2 )(2n+l + (? (A-^ 



f 


1/2- 







(17) 



This expression recovers the corresponding formula for the Hermitian generalised 
Ginocchio potential for q = 1. 
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Similarly to the case of the Hermitian version of the generalised Ginocchio potential 
the energy eigenvalues are written as E nq = — 7 4 /i^ g , and they are independent of e. 
Actually, we find that for the q = 1 choice the expressions for the Hermitian problem 
are recovered formally. However, the forthcoming analysis will show that despite the 
similar form, some quantities can be chosen complex for the VT symmetric case. Before 
going on we may note that fj, nq , and consequently E nq depends on the 2n + 1 + q(X — ~) 
combination, and this leads to a degeneracy between levels with q = 1 and q — — 1 
whenever A is a real half-integer number. Furthermore, for A = | states with opposite 
quasi-parity and with the same n become degenerate: in fact, this is the point where 
the spontaneous breakdown of VT symmetry sets in if the A parameter is continued to 
complex values allowed by Q. 

Let us now analyse the conditions for having real and complex energy eigenvalues 
E nq , which corresponds to inspecting the nature of /i rag fll7j) in terms of the allowed 
values of s and A displayed in Q. These also have to be combined with the condition 
Re(/i n q) > which guarantees normalisability of the solutions (fTH|) . The key element 
of the analysis is the term containing the square root in (fTTj) . so it is useful to inspect 
separately the cases when it is real, imaginary or complex, which corresponds to A > 0, 
A < and complex A, where 

A = 7 2 (, + \f + (1 - 7 2 ) (2n + 1 + q(X - \)j . (18) 

We restrict our analysis to 7 2 > 1: the alternative choice, 7 2 < 1 would change the 
nature of the r(u) function in (J2J). We can note that A occurs everywhere in the 
combination q(X — |), so when A is real, we can assume that A > |, because the A < | 
cases can be obtained simply by switching q = +1 to q = —1. Also, when A = | + il, it 
is enough to assume I > for the same reason. 

• A > 0. This can happen only if A is real, while from A > in (JTSj) and 7 2 > 1 
it follows that s also has to be real. Inspecting the allowed values of n for various 
parameter domains we find the following. Normalizable states can be obtained for 
\L nq in (JT7J) when 

-\ (l + 9(A - \) + (^) V2 1. + jl) < n < -i (l + 9 (A - \) - \s + i|) (19) 

holds. If the upper boundary of this domain is negative, then there are no 
normalizable solutions. This depends on the relative magnitude of s and A. 

• A < 0. Here again A has to be real, while s can take both real and complex values 
allowed in (jUJ). The Re(/x ng ) > condition now reduces to 2n + 1 + q(X — ~) < 0, 
which has to be combined with A < 0. The resulting condition is then 

n<-i(l +9 (A-i) + (^)"Vi|) (20) 

for real values of s and 

n^-^l + qiX- 1 -)) (21) 
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for s = — ~ + i<7. Note that these conditions can be met only for q = —1 if A is large 
enough (and positive, as we assumed before). 

• A is complex. For this A = \+il is required, while s can be both real and s = — ^+icr. 
This situation corresponds to the spontaneous breakdown of VT symmetry, and 
the energy eigenvalues appear in complex conjugated pairs due to (fJ> nq )* = ^ n ~qi 
which leads to {E nq )* = E n _ q . At the same time the Ke(fj, nq ) > condition turns 
out to be the same for q — +1 and —1, in accordance with the expectation that 
the number of normalisable states has to be the same for both quasi-parities. The 
detailed analysis is more complicated for complex values of A (and A) than for real 
A, so we can resort only to numerical calculations in this respect. The outcome 
depends on the relative magnitude of \s + \ \ (which is |er| for complex values of s) 
and /. 

We can now address the question whether with the tuning of the potential 
parameters the spontaneous breakdown of VT symmetry (i.e. the appearance of 
complex conjugate pairs of eigenvalues) happens at the same time for all the bound 
states as in the case of shape-invariant potentials [201 , or gradually, as for some non- 
shape-invariant potentials, such as the VT symmetric square well [23- From the analysis 
above we find that when this mechanism is realised via setting A to the complex value 
A = | + il, then all the energy eigenvalues turn to complex at the same time. The 
spectrum can also be changed to complex by tuning s from real to imaginary values 
and keeping A real. Since s is contained in the formulae in the combination s(s + 1) 
or (s + |) 2 = s(s + 1) + | which is always real, this possibility is more limited. In 
this case again all the energy eigenvalues turn into complex at the same time, but the 
character of the potential also changes, as its leading term in (0) changes sign. In the 
Hermitian setting this would correspond to replacing the potential well with a barrier, 
which obviously changes the nature of the problem. Although in the VT symmetric 
version of (0) the situation is less transparent, the complexification of the spectrum via 
tuning s to complex values and keeping A real is clearly different from the situation when 
s is kept real and the spontaneous breakdown of VT symmetry is induced by tuning A 
to complex values. 

In Figures 121 El and H] the real and imaginary components of (^) are plotted for 
fixed values of e, 7 and s and for various values of A corresponding to unbroken and 
spontaneously broken VT symmetry. The position of the energy eigenvalues are also 
indicated. 

A similar analysis can be performed for 7 2 < 1 too. In this case the s = — \ + ia 
choice plays a more important role, but otherwise the results are qualitatively the same. 
Note that in this case r(u) in (J2J) changes, and this also modifies the nature of the 
potential. 

Before closing this section we mention briefly some aspects of the one-dimesional 
version of the Ginocchio potential [2H] which is obtained from (Q) by the A = or 1 
substitution. This limit is analogous to the one-dimensional version of the harmonic 
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oscillator, which is also obtained from the radial harmonic oscillator after cancelling 
the singular centrifugal term, allowing the extension of the potential to the full x axis. 
Another similarity between the two systems is that the two choices of A correspond to the 
even and odd solutions, and this can clearly be seen from the structure of the bound-state 
solutions (jOJ, in which the Jacobi polynomial reduces to an even and odd Gegenbauer 
polynomial for A = and 1, respectively [37] . Losing the A parameter means that in 
the case of the VT symmetric one-dimensional Ginocchio potential the spontaneous 
breakdown of VT symmetry cannot be implemented as in the general case. It is also 
interesting to note that the quasi-parity quantum number occurs only in the combination 
q(\ — ~), which means that the q — +1, A = combination is equivalent with q = — 1, 
A = 1, and q — +1, A = 1 is equivalent with q = — 1, A = 0, and this reflects the relation 
of the quasi-parity quantum number with ordinary parity, similarly to the case of the 
one-dimensional harmonic oscillator [THj . It is also worthwhile to note that the Hermitian 
version of the one-dimensional Ginocchio potential possesses a number of complex- 
energy solutions (resonances) [26J, and since the energy eigenvalues are not sensitive 
to the e parameter appearing in the complex coordinate shift in (|SJ), these remain 
unchanged after the VT symmetrisation of the potential. However, these are unbound 
solutions, so their character is different from that of the (normalisable) complex-energy 
solutions which appear when the VT symmetry is spontaneously broken. 



4. Supersymmetric aspects of the VT symmetric generalised Ginocchio 
potential 



According to Ref. [23] the supersymmetric partner of a VT symmetric potential depends 
on the quasi-parity q, and thus corresponds to two distinct potentials. For this, the 
partner potentials have to be constructed by using a (/-dependent factorization energy 
such that |2S1 

Vi q \x) = U ( £{x) + e (q) = [W iq \x)] 2 ± + e (9) , (22) 

where = E$_ is the ground-state energy of the 'bosonic' potential , which, due 
to this construction is independent from q, i.e. Vl (%) = V(x). The superpotential 
W <yq \x) is expressed in terms of the ground-state (n = 0) wavefunction of V(x) ((TJ 

W i - q \x) = -^ln^(^) 

= 7 2 ( 7 2 - l)sinlm + 7 2 /j 0g sinh u 7 2 (| + g(A - \)) 

2(7 2 + sinh 2 w) 3 / 2 (7 2 + sinh 2 w) 1 / 2 (7 2 + sinh 2 u) 1 / 2 sinh u ' 
and it clearly depends on q explicitly (in third term) and implicitly via fj,Q q (in the 
second term). The 'fermionic' partner potentials of the generalised Ginocchio potential 
contain the same terms as (JIJ, but the coupling coefficients are different, and pick up 
g-dependence, as expected |2~H] : 

V to\ x) = ds + Bq coth2u + ^ W-i) 2 (24) 

7 2 + sinh 2 u 7 2 + sinh 2 « (7 2 + sinh 2 u) 2 4(7 2 + sinh 2 u) 3 ' 
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(25) 
(26) 



(27) 



These partner potentials are VT symmetric if fio q (and thus A too) are real, which, 
under rather general conditions, coincides with the requirement of the (unbroken) VT 
symmetry of ([IJ itself. This was the case for some VT symmetric shape-invariant 
potentials too |2lH 12*2*] . When A = ~ + il, which happens when the VT symmetry of 
is spontaneously broken, (J2"3j) ceases to be VT symmetric, which is again a result 
similar to those obtained for shape-invariant potentials Pffi I22j. 

Figures El El and E| display also the 'fermionic' partners of the respective 'bosonic' 
potentials. Due to the SUSYQM construction the energy eigenvalues of the 'fermionic' 
partners are the same with the exception that the levels with n = and q — ±1 are 
missing from the spectrum of V± (x) . The example in figure |U corresponds to the 
spontaneous breakdown of the VT symmetry of the 'bosonic' potential, and thus the 
VT symmetry of the 'fermionic' potentials is manifestly broken. This is indicated by the 
fact that the real and imaginary component of the potential ceases to have definite parity 
under space reflexion. However, the two 'fermionic' potentials are the VT transforms 
of each other: V+ (x) = [vj. ^(— x)]*. Note that for A = |, i.e. for the point of the 
spontaneous breakdown of VT symmetry the two 'fermionic' partners with q = +1 and 
q = —1 coincide. 

5. Summary and conlcusions 

We analysed a Natanzon-class potential, the generalised Ginocchio potential in a VT 
symmetric setting in order to explore similarities and differences with the more restricted 
shape-invariant potential class. This work was inspired by the fact that up to now 
the exactly solvable VT symmetric potentials were almost exclusively members of 
the shape-invariant class, and they showed marked differences compared to examples 
outside this class, e.g. those which have been solved numerically or had quasi-exactly 
solvable character. Our analysis showed that the VT symmetric generalised Ginocchio 
potential shares all the specific properties of shape-invariant potentials. In particular, 
its states can also be characterised by the quasi-parity quantum number, and the 
spontaneous breakdown of its VT symmetry takes place suddenly, i.e. by tuning a 
potential parameter (A) all its real energy eigenvalues turn into complex conjugate 
pairs at the same value of this parameter. These results seem to originate from the 
'robust' structure of the normalisable solutions of Natanzon-class potentials, which 
allows the implementation of VT symmetry to these technically non-trivial problems. 
Another similarity with the shape-invariant potentials is that the VT symmetric 
generalised Ginocchio potential has two 'fermionic' supersymmetric partners (generated 
by eliminating the lowest state of the original 'bosonic' potential with quasi-parity 
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q = +1 and —1), and the partner potentials also possess VT symmetry as long as 
the VT symmetry of the 'bosonic' potential is unbroken, but they cease to be VT 
symmetric when the VT symmetry of the 'bosonic' potential is spontaneously broken. 
This seems to indicate that the properties thought to be specific to VT symmetric 
shape-invariant potentials might be valid to the much larger Natanzon potential class 
too, and perhaps also beyond that. Further studies should be made to check the validity 
of this conjecture. 
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Figure 1. The real (left panel) and imaginary (right panel) component of the u(x) 
function for 7 = 1.75 and e = 0.3. Note the different vertical scales. 




Figure 2. The real (left panel) and imaginary (right panel) component of the potential 
||TJ for e — 0.3, 7 = 1.75, s = 8.1 and A = 1.25 (solid line) and its supersymmetric 
partners Vj +1 \x) (dashed line) and vj_ ^(x) (dotted line) in l)24[l. Normalisable states 
of (HJ arc found at E +1 = -171.313, E x +1 = -106.160, E 2 +i = -46.679, E 3 +1 = 
-5.666; E _i = -218.913,^1 _ x = -154.978, E 2 -1 = -90.379, E 3 _ x = -33.993 
and E4 _i = —1.061. The spectrum of V±(x) is the same, with the exception of the 
Eq q level, which is missing from its spectrum. 
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Figure 3. The same as Figure [21 with A = 0.5. Normalisable states are found at 
E +i = E _i = -195.477, E x +1 = E x _ x = -130.419, £ 2 +x = £ 2 -1 = -67.675 
and £3 = E 3 _i = —17.640. The supersymmetric partners vl +1 \x) and V+ X \x) 
coincide in this case. 




Figure 4. The same as Figure|3]with A = 0.5 + 1.25i corresponding to spontaneously 
broken VT symmetry. Normalisable states are found at Eq +1 = (Eq = 
-196.494 + i 40.038, E x +1 = (E x _ x )* = -130.023 + i 41.130, £2 +1 = (E 2 -1)* = 
-65.367 + i 37.105 and E 3 +1 = (E 3 = -11.833 + i 25.161. The supersymmetric 
partners Vj_ (x) and vl ^(x) cease to be VT symmetric in this case. 



